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LASER  DRIVEN  ELECTRON  ACCELERATION  IN 
VACUUM,  GASES  AND  PLASMAS 


I.  INTRODUCTION 

Breakdown  and  thermal  limitations  in  the  cavities  of  radio  frequency 
linacs  limit  the  accelerating  gradient  to  <  50  MV/m.  Intense  focused 

12 

optical  beams,  on  the  other  hand,  can  reach  field  levels  greater  than  10 

V/m  in  vacuum.  Such  enormous  fields  have  stimulated  a  great  deal  of 

research  in  laser  driven  acceleration  concepts.^  This  paper  will  discuss 

2-8 

laser  acceleration  of  electrons  in  i)  vacuum,  ii)  neutral  gases  (inverse 

9-11 

Cherenkov  acceleration),  and  iii)  plasmas  (laser  wakefield 

,  V  12-25 

acceleration). 

26  27 

Advances  in  laser  technology  ^  have  made  possible  compact  teravatt 

18  2 

laser  systems  with  high  intensities  (>  10  V/cm  ),  modest  energies  (<  100 
J),  and  short  pulses  (<  1  psec).  The  peak  amplitude  of  the  transverse 
electric  field  of  a  linearly  polarized  laser  pulse  is  given  by 

E,  [TV/m]  =  2.7  x  10"^I^^^[W/cm^]  =  3.2  a  /X[vm],  (1) 

where  I  is  the  laser  intensity,  X  is  the  laser  wavelength,  and  a^  is  the 

laser  strength  parameter.  Physically,  a^  is  the  normalized  transverse 

oscillation  momentum  of  the  electrons  in  the  laser  field.  The  power  of  a 

2 

linearly  polarized,  Gaussian  laser  pulse  is  P[GW]  =  21.5(a^r^/X)  ,  where  r^ 
is  the  minimum  spot  size. 


II.  LASER  ACCELERATION  IN  VACUUM 


2-8 

The  acceleration  of  electrons  in  vacuum  by  optical  fields  is 

limited  by  diffraction  effects  and  by  electron  slippage.  The  phase 

velocity  of  the  optical  field  co-propagating  in  the  direction  of  the 

accelerated  electrons  is  greater  than  c  and  given  approximately  by  v 

2 

1  +  l/(kZj^),  where  «  =  ck  is  the  laser  frequency,  =  kr^/2  is  the 


ph 


/c  = 
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Rayleigh  length,  i.e..  Z,  la  the  distanee  over  .hich  the  spot  sire  expands. 
Since  V  .  >  c.  electrons  vith  =  c  will  phase  slip  with  respect  to  the 
accelerating  field  and  decelerate.  This  will  occnr  over  a  slippage 
distance  z^,  which  lor  highly  relativistic  electrons  is  -  Z,(. 

A  simile  calculation  based  on  Maxwell's  equations  in  vacuum  shows  that 
under  certain  restrictive  conditions  no  net  energy  gain  is  possible  using 
optical  fields.  This  has  become  known  as  the  Lawson-Woodward  <L») 
theorem.2-8  ^he  LV  theorem  assumes:^'^  (1)  the  region  of  interaction  is 
infinite,  (li)  the  laser  fields  are  in  vacuum  with  no  walls  or  boundaries 
present,  (iii)  the  electron  is  highly  relativistic  (v  =  c)  along  the 
acceleration  path,  (Iv)  no  static  electric  or  magnetic  fields  are  present, 
and  <v)  nonlinear  effects  (e.g.,  ponderomotive,  v  x  B,  and  radiation 
reaction  forces)  are  neglected.  One  or  more  of  the  assumptions  of  LU 
theorem  must  be  violated  in  order  to  achieve  a  nonzero  net  energy  gain. 

For  example,  finite  energy  gain  can  be  achieved  by  introducing  a  background 
of  gas  into  the  interaction  region,  as  in  the  inverse  Cherenkov 
accelerator.’-^^  The  gas  can  reduce  the  phase  velocity  of  the  optical  beam 
to  less  than  c,  reducing  the  slippage,  but  diffraction  remains  a 
limitation.  Alternatively,  acceleration  can  result  fro.  the  introduction 
of  a  periodic  magnetic  wiggler  field,  as  in  the  inverse  free  electron 
laser. 2®-^*  In  vacuum,  a  nonzero  energy  gain  can  be  achieved  by  the 
introduction  of  boundaries  which  limit  the  interaction  distance  to  a  region 

r~  z  >  about  the  laser  focus. 

'  s 


A.  ArrPleration  Using  Higher  Order  Gaussian  Beajns 

Higher  order  Gaussian  inodes  can,  in  principle,  provide  an  axial 
electric  field  component  for  electron  acceleration  in  vacuum.^  By 
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properly  choosing  the  electron  injection  point,  net  acceleration  is 
possible.  Consider  a  radially  polarized  higher  order  Gaussian  mode 
propagating  along  the  z-axis  in  the  positive  z-direction,  after  having  been 
reflected  off  a  mirror  located  at  some  negative  z-position,  as  shown  in 


Fig.  1.  The  results  can  be  readily  generalized  to  describe  more 

complicated  configurations.  Near  the  axis,  E^.  =  Ej^(rrp/rg)sin<j/  and  E^ 

22  1/2  . 

E, (2r  /kr^)cos\|>,  where  =  kz  -  «t  -  2tan  r^  =  r^(l  +  z  /Zj^)  is 

1  o  s 

the  laser  spot  size,  =  nr^/X,  and  E^  is  a  constant.  The  phase  velocity 
along  the  axis  (r  =  0)  and  near  the  focal  point,  |z|  <  Zj^,  is  given  by 

ph' 


v_  /c  =  1  +  1/(2y^),  where  y  =  nr  /42X  defines  a  critical  energy.  The 


pil  V.  w 

slippage  distance,^’®  z  ,  defined  as  the  distance  required  for  the  electron 


to  phase  slip  by  n,  is  given  by  wz^ |v 


-1  -1 
-  V 


1  =  n. 


Zg  =  (nZj^/2)(l  +  y^/y^) 


(2) 


where  y  =  (1  -  v^/c^)"^'^^.  In  the  high  energy  limit  (y  »  Zg  =  nZj^/2 

and  in  the  low  energy  limit  (y  «  Y^,),  z^  =  X^  «  Z^.  Figure  2  shows  the 
accelerating  field  as  a  function  of  axial  distance  for  high  (solid  curve) 
and  low  (dashed  curve)  energy  injection.  In  both  cases,  the  total  area 
under  the  curve  is  zero,  i.e.,  there  is  no  net  energy  gain  from  -®  <  z  <  ». 

For  a  finite  interaction  region  centered  about  the  origin,  the  maximum 
energy  gain^’®  AW  is  given  approximately  by  the  peak  amplitude  of  the  axial 
electric  field,  2E^/kr^,  multiplied  by  the  slippage  distance,  z^,  i.e., 

AWtMeV]  =  31P^^^[TW](1  +  y^/y^)~^  <3) 

2 

where  P  =  cE^r^/32  is  the  laser  power.  The  critical  energy  =  me  y^,»  cari 
be  written  as  W^[MeV]  =  l.l(r^/X).  In  the  high  energy  limit  (y  »  Y^),  the 
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energy  gain  can  be  substantial,  i.e.,  AW  =  100  MeV  for  P  =  10  TW.  In  the^ 
low  energy  limit,  however,  this  energy  gain  is  reduced  by  the  factor  y  /y^ 
«  1.  When  damage  thresholds  are  considered,  the  low  energy  limit  y  «  Y^, 
appears  to  be  the  relevant  regime  for  typical  parameters  of  interest. 

B,  Limitations  Due  to  Material  Damage 

In  principle,  limiting  the  interaction  distance  to  a  small  region  near 

the  focus  can  lead  to  substantial  energy  gains.  In  practice,  however,  the 

energy  gain  can  be  limited  by  the  intensity  damage  threshold  of  the 

reflecting  surface  material.®  As  an  example,  consider  placing  a  mirror  at 

a  distance  -z  /2  from  the  focus  (z  =  0)  and  using  the  higher  order  Gaussian 
s 

mode,  see  Fig.  1.  The  electron  energy  gain  is  one-half  the  value  given  by 
Eq.  (3). 

The  laser  intensity  on  the  surface  of  the  mirror  1^  must  be  less  than 

the  damage  threshold  limit,  I^  =  P/dtr^)  <  I^,  where  r^  is  the  radiation 

spot  size  on  the  mirror  surface  and  I^  is  the  mirror  damage  threshold 

intensity.  Typically, for  a  1  ps  laser  pulse,  1^  <  ^  TW/cm^.  Since 

z  <  Z  ,  r  =  r  .  Hence,  for  a  fixed  laser  power,  1^  <  I^j  implies  r^  > 
s  ^  R  m  o 

(P/nI^)^^^.  Increasing  r^  increases  the  critical  energy,  since  y^, 
nr  /(^X).  The  condition  1^  <  1^  implies  y^  >  iiP/(2I^X^),  which 
coLsponds  to  a  critical  energy  W^  =  mc^^,  W^IGeV]  =  (6.4/Xlum])(PlTWl/ 
I^[TW/cm2])^^^.  Typically,  this  value  of  W^  is  quite  high.  For  P  =  10  TW, 
X  =  1  ym,  and  Ij  =  5  TW/cm^,  the  injected  beam  energy  should  be  greater 
than  W  =  9  GeV  to  be  in  the  high  energy  limit.  On  the  other  hand,  if  the 
injected  energy  is  below  this  critical  value,  the  energy  gain  is  one-half 
that  given  by  Eq.  (3)  with  y  «  Y^-  For  X  =  1  ym,  P  =  10  TW,  1^  =  5  TW/cm 
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and  Wj  =  1  GeV,  the  energy  gain  is  small,  AW  <  0.6  MeV.  These  same 
arguments  regarding  the  damage  threshold  can  be  applied  to  other  laser 
acceleration  configurations  in  vacuum  which  rely  on  limiting  the 
interaction  region  by  using  optical  components. 


C.  Limitations  Due  to  Apertures 

To  limit  the  electron  and  optical  beam  interaction  distance  within  an 
acceleration  stage,  the  electron  bunch  must  propagate  through  an  aperture 
and  the  optical  beam  reflected,  see  Fig.  1.  A  portion  of  the  self-fields 
associated  with  the  electron  bunch  extends  beyond  the  aperture  and  will 
essentially  be  reflected  by  the  boundary  (dielectric  mirror).  A  limit  on 
the  maximum  accelerated  charge  per  bunch  is  reached  when  the  reflected 
self-field  energy  equals  the  energy  gain  per  stage  or  when  the  self-fields 
on  the  reflecting  boundary  become  comparable  to  the  damage  threshold  value 
The  self-fields  associated  with  a  highly  relativistic  electron  are 


dominated  by  the  radial  electric  field  E^.  and  azimuthal  magnetic  field  B^, 

i.e.,  the  radial  electric  field  is  larger  than  the  axial  electric  field  by 

the  factor  y^-  The  self-fields  for  a  single  electron  of  charge  q,  at 

position  r  =  0,  and  2  .  .  v^t.  Is  -  qrrlr  .  t  (s-z„)  )  mi 

R  -  (V  /Cl  E  , .  The  radial  field  of  an  electron  bunch  of  length 
“0,1  ■  ''  r,l 

radius  r.  and  consisting  of  N  electrons,  evaluated  at  z  =  0  (midplane  of 
D 

bunch)  and  r  >  is 


2o2,,^-l/2 


h  2  2  2 


(4) 


-t./2 


The  total  self-field  energy  which  intercepts  the  boundary  is  given 
approximately  by 
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00  2j^/2 

"self  “  h  J 
"a 

\ 

=  (qV/£b)[2tn((l  +  -  tan'^  ,  (5) 

where  Zj,  =  ^b  aperture  radius,  r^  »  rj^,  and  = 

2r  /(Y^i,).  The  fact  that  the  radial  self-field  of  a  single  electron  is 
confined  within  an  angle  1/y  is  apparent  in  the  limits  of  integration  over 
2  in  Eq  (5).  The  self-field  energy  intercepting  the  aperture  in  the 
extended  electron  bunch  limit,  »  Zr^/Yf  is 

'^self  =  (‘^^N2/e^^)(2en(YV2»^a> 

and  in  the  localized  electron  bunch  limit,  «  2r^/Y>  is  '^ggif 
qVY/2r^. 

Equating  the  total  energy  gain  per  stage  (NAW)  with  the  total 
reflected  self-field  energy  (RW^^if)  places  a  limit  on  the  amount  of  charge 
which  can  be  accelerated,  where  tV  is  the  single  electron  energy  gain  per 
stage  and  R  is  the  self-field  energy  reflection  coefficient.  The  ^ 

reflection  coefficient  is  approximately  given  by  R  =  jd  -  ’re)/(l  +  ’Te)  1  . 
where  e  is  the  dielectric  constant  of  the  reflecting  surface  and  is  taken 
to  be  independent  of  frequency  even  though  the  self-fields  have  a  frequency 
spectrum  which  peaks  around  =  nc/£j,.  In  the  extended  electron  bunch  limit 
the  maximum  number  of  electrons  which  can  be  accelerated  to  energy  Ymc^  is 

.-1 

N  =  (iu/r  )(AY/2R)(en(e.  Y/2r  )  +  hMJ  ,  (7) 

max  'be  v  u  a 
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and  in  the  localized  electron  bunch  limit  is  =  (2r^/rg)( Ay/y)/R  where 

2 

r  =  a^/mc^  is  the  classical  electron  radius  and  Ay  =  AW/mc  .  For 
e  ^ 

highlyrelativistic  beams,  the  extended  bunch  limit  is  relevant,  >> 

4 

As  an  example,  consider  a  10  GeV  electron  bunch  (y  =  2  x  10  )  with  =  0.1 
wm.  Ay  =  1,  R  =  0.5,  and  r^  =  5  ym.  Equation  (7)  yields  =  6  x  10^, 

whereas  the  peak  self-field  at  the  edge  of  the  aperture  is  =  30  GV/m,  as 
obtained  from  Eq.  (4)  with  N  =  6  x  10^.  This  self-field  level  exceeds 
typical  damage  threshold  values,  which  are  ~  1  GV/m. 

D.  Vacuum  Beat  Wave  Accelerator 

The  vacuum  beat  wave  accelerator  (VBWA)  is  an  example  of  an 

7  8 

accelerator  based  on  the  nonlinear  ponderomotive  force.  ’  In  the  VBWA, 

two  laser  beams  of  different  wavelengths,  and  X2,  are  co-propagated  in 

the  presence  of  an  injected  electron  beam.  Properly  phased  electrons, 

traveling  essentially  along  the  same  axis  as  the  two  laser  beams, 

experience  an  axial  acceleration  from  the  beat  term  in  the  v  x  B  force.  By 

adjusting  the  wavelengths  and  spot  sizes  of  the  laser  beams,  the  phase 

velocity  of  the  ponderomotive  beat  wave  can  be  adjusted  such  that  v^j^  <  c. 

The  acceleration  mechanism  in  the  VBWA  is  similar  to  that  of  the  inverse 

7  8 

free  electron  laser. The  energy  gain  in  the  VBWA  is  given  by 

Wp[MeV]  =  [Wj[MeV]  +  480(Xj/X2  -  l)Pj[TW]]^^^.  (8) 

where  W_  (W-)  is  the  final  (initial)  electron  energy  and  we  assumed  Wj  » 

F  j. 

mc^.  As  an  illustration,  for  Xj^  =  2X2  =  1  um  and  Pj^  =  20  TW,  the  energy 
gain  is  100  MeV. 
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III.  laser  acceleration  in  a  gas 

The  linear  dispersion  relation  for  a  laser  beam  propagating  in  a 
neutral  gas  implies  an  axial  phase  velocity  given  by  Vpj^/c  -  1  + 

_  i)/2,  where  is  the  linear  refractive  index.  Typically,  -  1  « 
1  and  is  proportional  to  the  neutral  gas  density  n^^.  Proper  choice  of  n^^ 
can  result  in  <  c  and  electron  phase  slippage  can  be  eliminated, 
however,  diffraction  remains  an  important  limitation.  This  is  the  basis  of 
the  inverse  Cherenkov  accelerator  (ICA),’-“  vhich  uses  a  higher  order 
Gaussian  laser  beam  to  accelerate  electrons  in  a  gas.  Experiments  at  BNL 
on  the  conventional  ICA  observed  a  3.7  MeV  energy  gain  (31  MeV/m)  of  an 
injected  electron  beam  <40  HeV)  using  a  580  HW  COj  laser  in  2.2  at.  of  Hj 

gas. 


A.  Self-Guided  Inverse  Cherenkov  Acceleration 

In  the  conventional  ICA,  the  electron  accelerating  distance  (energy 
gain)  is  severely  limited  by  diffraction  effects.  By  self-guiding  the 
optical  beam  in  the  ICA,^^  substantially  higher  electron  energy  gams  can 

be  achieved. 

A  self-guided  inverse  Cherenkov  accelerator  (ICA)  has  recently  been 

proposed  and  analysed.*^  The  self-guided  ICA  uses  a  higher  order  Gaussian 

laser  pulse  vhich  has  a  radially-polarised  electric  field  component  of  the 

form  E  m  E  <42r/r  )exp(-r2/r2  .  io)  .here  r^  is  the  laser  spot  sise  and  ♦ 
r  o  s  ^ 

is  the  laser  phase  (near  the  axis,  ♦  .  ks  -  mt  -  2tan  (z/Tp))-  1" 

addition,  there  is  a  significant  accelerating  axial  field  peaked  along  the 
z-axis  vith  magnitude  \eJ  =  2j2E„/kr^  .here  k  -  2n/X.  The  self-guided  ICA 
configuration  operates  at  laser  po.ers  near  the  critical  peer  for  self- 
focusing  and  at  intensities  for  vhich  the  gas  is  slightly  ionised.  Figure 
3  shove  the  radial  profiles  of  the  electric  fields  and  plasma  for  the  self- 
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guided  ICA. 

The  propagation  of  intense  optical  beams  in  gases  is  affected  by  a 
combination  of  diffraction,  refraction,  and  ionization.  The  refractive 

0  9  y/  2  /  vl/2 

index  Is  given  by  n(r)  =  *  -IjKD  -  (o;(r)/2»  ,  where  .  (4llq  n^/m) 

is  the  electron  plasma  frequency,  n^  is  the  electron  plasma  density,  ^2  is 
the  nonlinear  refractive  index,  and  I  is  the  intensity  of  the  optical  beam. 
Generally,  ^2  is  positive  and  results  in  self-focusing  of  the  optical  beam 
when  the  beam  power  is  greater  than  the  nonlinear  focusing  (critical) 
power,  which  for  a  fundamental  Gaussian  beam  is  =  X  /(2iiri^»l2) •  As 
the  beam  self-focuses  the  peak  intensity  becomes  sufficiently  large  so  that 
tunneling  ionization  occurs  and  a  plasma  is  generated.  The  local  decrease 
in  the  refractive  index  due  to  the  plasma  tends  to  defocus  the  optical 
beam.  If  diffraction,  self-focusing  due  to  h2  defocusing  due  to. plasma 
generation  are  properly  balanced,  a  self-guided  optical  beam  can  be  formed. 

An  envelope  equation  describing  the  evolution  of  the  optical  beam  spot 
size  r  (z,t)  can  be  derived  by  applying  the  source  dependent  expansion 
(SDE)  method^^  to  the  paraxial  wave  equation.  In  the  SDE  method,  the 
optical  beam  is  expanded  in  a  complete  set  of  orthogonal  Laguerre-Gaussian 
functions  which  are  implicitly  functions  of  C  =  z  -  ct  and  the  propagation 
time,  T  =  t,  through  the  optical  beam  parameters,  i.e.,  spot  size, 
wavefront  curvature,  amplitude  and  phase.  The  self-guiding  analysis 
assumes  that  the  optical  beam  is  adequately  described  by  a  single  Laguerre- 
Gaussian  node  and  higher  order  node  coupling  and  generation  are  neglected. 
Ionisation  is  considered  in  the  high  field  limit  and  is  nodeled  by  tbe 
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tunneling  ionization  rate.^^  The  evolution  of  the  spot  size  rg(C,T)  is 
w  11 

given  by 

32r/3,2  .  (c/Zj„)V3(l  -  P/P^  ^  <’> 

where  R  =  r  /r  is  the  spot  size  normalized  to  its  initial  value,  = 

kr2/2  is  thi  iLtial  Rayleigh  length,  P  -  (cnyiSiE^rJ  in  the  total  peer, 
o 

p  =  8n/(k.^Vl  ho)  is  the  nonlinear  focusing  power  for  the  radially 
polarized  beam,  =  Anq^/mc^,  and  is  the  filling  factor.  The  terms 
on  the  right-hand  side  of  Eq.  (9)  denote  respectively,  vacuum  diffraction, 
nonlinear  focusing  and  plasma  defocusing.  The  filling  factor  is  given  by 


=  J  dx(n^/nj^)(2  -  x)xe 

o 


(10) 


where  x  =  Zr^/r^.  In  the  limit  n  /n^  «  1,  the  plasma  density  is  given  by 
s 

an^/aC  =  -  (V^^'^ion^  '^ion^l-l^  tunneling 

ionization  rate,^^  which  is  a  function  of  the  field  amplitude  jEl- 

The  condition  for  a  matched  beam,  i.e.,  a^R/ax  =  0,  is  given  by 

1  ^  1^2  j.2^  ^  Yor  a  matched  beam,  it  can  be  shown  that  the  phase 

pn  o  R  -  2 

velocity  is  given  by  Vpj^/c  =  1  -  0.5(X/nr^)  -  (*^0  "  l)/2.  Note  that 

V  <  c  and  is  determined  by  the  linear  and  nonlinear  refractive  indices  as 
ph 

well  as  the  power. 

The  accelerating  gradient  |Ej  -  (^/n)(X/rpE„  can  be  estimated  by 
considering  the  case  of  a  matched  beam  i.e.,  P  =  P„[_  and  >  0.  For  this 
case  ve  find  that  E„  »  (8nel/c)l'^  Vr^  =  (n/aXeb^hjI)^^^  and  the 
accelerating  gradient  becomes 


|E  I  .  (eln^ll2^/2)^'\  =  2e(llt,^n2ll/c)‘'‘Ii 


1/2.,  1/2 


(11) 
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where  we  have  set  I  equal  to  the  ionization  intensity  (I  =  Ij)- 

The  phase  velocity  can  be  controlled  by  introducing  a  small  amount  of 
background  plasma. A  transversely  uniform  background  plasma  will 
increase  the  phase  velocity  but  have  no  effect  on  the  focusing  of  the 
optical  beam.  The  background  plasma  can  be  created  by  introducing  a  small 
concentration  of  easily  ionized  gas,  i.e.,  low  ionization  intensity.  In 
addition,  the  background  density  can  be  tapered  as  a  function  of  z  to 
increase  the  phase  velocity  and  optimize  electron  acceleration. 


B.  lonization-Modulation  Instability 

The  equilibrium  described  above  is  subject  to  an  ionization-modulation 
(IM)  instability.^^  The  IM  instability  is  due  to  varying  degrees  of 
ionization  along  the  beam  and  results  in  the  modulation  of  the  optical  beam 
envelope.  The  IM  instability  can  be  analyzed  by  perturbing  the  equilibrium 
spot  size  of  an  optical  beam,  R  =  +  SR,  where  R^  is  the  equilibrium 

value  and  SR  is  the  perturbation.  Expansion  of  Eq.  (9)  about  R^  indicates 
that  the  perturbed  envelope  satisfies  an  equation  of  the  form 

.  J°u(c')®(i;’)<ic'.  (12) 

c 

where  U  is  a  slowly  varying  function  of  C*  For  a  fundamental  Gaussian 
beam,  it  can  be  shown^^  that  in  the  linear  regime  the  asymptotic  growth  of 
the  IM  instability  is  of  the  form  SR  -  exp[(l  ±  i’r3)Ng],  where 

1/3 

The  growth  of  the  IM  instability  depends  on  the  profile  of  the  plasma 
electron  density  along  the  axis,  As  an  example,  for  a  gas  at  20  atm 
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with  n  /n  =  10’^  and  r  «  28  wm,  as  in  Fig.  4  discussed  below,  the  beam 
undergoes  approximately  25  e- folds  (N^  ~  25)  after  propagating  a  distance 

of  100  Rayleigh  lengths  (ex  =  100  Z^^). 

To  gain  some  understanding  of  the  IM  instability,  consider  increasing 
the  spot  size  of  a  matched  optical  beam,  i.e.,  8R(x  =  0)  >  0.  In  this  case 
the  beam  intensity  and  ionization  rate  are  reduced  resulting  in  less  plasma 
generation  and  focusing  of  the  beam.  The  focusing  optical  beam  overshoots 
its  equilibrium  value  such  that  SR  <  0  some  distance  behind  the  beam  front. 
When  SR  <  0,  the  intensity,  ionization  rate,  and  plasma  density  increase, 
causing  the  beam  to  defocus  and  overshoot  its  equilibrium  value.  This 
focusing  and  defocusing  of  the  beam  results  in  the  IM  instability.  The 
modulation  amplitude  and  period  are  functions  of  the  distance  back  from  the 
head  of  the  optical  beam,  \l,\,  and  the  propagation  distance  cx. 

C.  Ionization  and  Collisional  Losses 

As  the  plasma  electrons  oscillate  in  the  laser  field,  they  can  collide 

with  the  background  plasma  electrons,  ions,  and  the  neutral  atoms.  For  a 

weakly  ionized  gas,  electron  scattering  is  dominated  by  electron-neutral 

collisions. The  electron-neutral  collision  frequency  is  given  by  = 

no  <v  >,  where  a  is  the  electron-neutral  cross-section  and  is  an 

n  en  os  '  en 

average  oscillation  velocity  of  the  electrons.  Typically,  is 

-15  -3 

approximately  constant  for  electron  energies  e  <  Eq,  where  ~  10  cm 

19  -3 

and  e  is  on  the  order  of  tens  of  eV.  For  example,  n^^  =  3  x  10  cm  and 

^  .  -1 
<v  >  =  5  X  10”^c  imply  an  electron-neutral  collision  time  of  x^^^  =  ~ 

os 

200  fs. 
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Ionization  of  the  gas  as  well  as  electron-neutral  collisions  result  in 
the  depletion  of  the  laser  intensity,  I,  which  is  found  to  be  given  by 
31/3(0-1)  =  where 


4.  U  k 

en  p  ion  pn 


(14) 


is  the  intensity  attenuation  length.  In  Eq.  (14),  is  the  tunneling 

ionization  rate,^^  kp  =  »p/c,  f  =  ‘^p  cross- 

sectional  area  of  the  generated  plasma  (laser)  and  the  oscillation  energy, 

in<v  >^/2,  is  assumed  large  compared  to  the  ionization  energy.  For  =  3 
OS  ^  « 

X  10^^  cm"^,  n^  =  10^^  cm"^,  =  5  x  10^^  s”^,  =  2x10  s  ,X=1 


pm,  and  f  ~  10  the  attenuation  length  is  ~  30 


m. 


D.  Numerical  Examples 

As  an  example,  consider  a  self-guided  ICA,  in  which  a  higher  order, 
radially  polarized  optical  beam  of  wavelength  X  =  1  ym  propagates  in 
hydrogen  (H2)  at  30  atm.  Since  E^  ~  r\Y^  -  n^^,  the  accelerating  gradient 
can  be  increased  by  increasing  the  gas  pressure.  Hydrogen  is  chosen  for 
its  low  atomic  number,  Z,  i.e.,  low  bremsstrahlung  losses  Z  )•  At  30 
atm,  n^^  =  8.1  x  10^°  cm"^,  ^2  “  3.3  x  10'^®  cm^/W,  =  1.9  GW  and 
ionization  potential  Uj.  =  15.4  eV.  The  matched  beam  conditions  can  be 
determined  from  Eqs.  (9)  and  (10)  for  a  given  axial  intensity  profile.  For 
a  constant  peak  intensity  of  I  =  4.7  x  10^^  W/cm^,  the  matched  profiles  for 

the  spot  size  r  and  the  power  P  are  shown  in  Fig.  4.  For  these 

®  -4 

parameters,  the  amount  of  ionized  plasma  is  small,  n^/n^^  <  10  .  The  on- 

axis  accelerating  field,  also  shown  in  Fig.  4,  has  the  maximum  value  450 

MV/m  at  the  front  of  the  optical  beam. 

As  discussed  above  the  optical  beam  undergoes  a  IM  instability. 

Numerical  simulations  of  Eq.  (9)  show  that,  with  a  1%  initial  perturbation 
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of  the  spot  size,  the  IM  instability  significantly  disrupts  the  beam 

envelope  after  ~  10  cm  of  propagation.  The  growth  rate  of  the  IM 

instability  is  a  highly  nonlinear  function  of  the  laser  intensity,  through 

13 

the  electron  density.  Reducing  the  intensity  in  the  example  to  3.2  x  10 
M/cm  (which  reduces  to  300  MeV/ra)  results  in  a  matched  beam  with  very 
little  ionization  (the  plasma  density  is  reduced  by  a  factor  of  130)  and 
little  variation  in  power  and  spot  radius  along  the  length  of  the  laser 
pulse.  Simulations  show  that  with  a  1%  perturbation  in  the  spot  size,  the 
pulse  propagates  >  1  meter  without  significant  disruption. 

A  self-focusing  mechanism  has  recently  been  reported  with  ultrashort 
fundamental  Gaussian  laser  pulses. In  these  experiments  laser  pulses 
with  intensities  of  ~  7  x  10^^  W/cm^  and  spot  sizes  of  -  40  ym  were 
propagated  over  20  m  in  air.  This  experiment,  however,  does  not 
demonstrate  the  self-guiding  of  a  fundamental  Gaussian  beam  since  a 
significant  portion  of  the  beam  energy  exists  beyond  the  40  ym  spot  size. 

IV.  LASER  ACCELERATION  IN  PLASMAS 

Plasmas  can  offer  some  advantages  as  an  accelerating  medium  in  laser 

driven  accelerators.  Plasmas  can  be  modulated  and  sustain  ultrahigh 
electric  fields,  and  under  appropriate  conditions  can  optically  guide  the 
laser  beam.  The  two  laser-plasma  accelerator  configurations  which  have 
received  the  most  attention  are  the  plasma  beat  wave  accelerator 

and  the  laser  vakefield  aeeelerator  (LWPA).!^-”  The  PEWA 
Utilizes  two  laser  beams  of  frequencies  and  such  that  -  «2  =  «p 
is  the  plasma  frequency,  to  resonantly  drive  a  plasma  wave.  The  following 

discussion  is  concerned  primarily  with  the  LWFA. 
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A. 


Laser  Wakefield  Acceleration 


As  an  intense  laser  pulse  propagates  through  an  underdense  plasma, 
i.e.,  X2/X2  «  1,  where  Xp  =  2nc/(0p,  (Op  =  (4nq\^/m)^^^  and  n^^  is  the 
ambient  electron  density,  the  ponderomotive  force  associated  with  the  laser 
pulse  envelope,  Fp  -  7a^  expels  electrons  from  the  region  of  the  laser 
pulse.  If  the  pulse  length  is  approximately  equal  to  the  plasma 
wavelength,  CT^  =  Xp,  the  ponderomotive  force  excites  large  amplitude 
plasma  waves  (wakefields)  with  phase  velocities  approximately  equal  to  the 
laser  pulse  group  velocity. The  maximum  wakefield  amplitude  generated 
by  a  linearly  polarized  laser  pulse  of  amplitude  a^,  in  the  1-D  limit  r^  » 
.2  .  15-18 

Xp,  IS 

EmaxlGeV/m]  =  3.8xlO"®(ng^Icm"^l )  a^(l  +  a^/2)  .  (15) 


In  the  absence  of  optical  guiding,  the  interaction  distance,  L.^^,  will  be 

limited  by  diffraction,  i.e.,  The  maximum  energy  gain  of  the 

electron  beam  in  a  single  stage  is  AW  =  which,  in  the  limit  aj  « 

1,  may  be  written  as^®  AW[MeV]  =  580(X/Xp)P[TWl .  As  an  example,  consider  a 

T  =  1  psec  linearly  polarized  laser  pulse  with  P  =  10  TW,  X  =  1  ym  and  r^ 
L 

=  30  ym  (a^  =  0.72).  The  requirement  that  ct^^  =  Xp  implies  n^^  =  1.2  x 
10^®  cm-3.  This  gives  =  2.0  GeV/m,  =  0.9  cm,  and  AW  =  18  MeV. 

The  interaction  length  and  the  electron  energy  gain  may  be  greatly 
increased  by  optically  guiding  the  laser  pulse  in  the  plasma. 


B.  Optical  Guiding  in  Plasma 

For  a  laser  pulse  to  propagate  in  a  plasma  beyond  the  limits  of  vacuum 
diffraction,  some  form  of  optical  guiding  is  necessary.  For  sufficiently 
powerful,  long  laser  pulses,  diffraction  can  be  overcome  by  relativistic 
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The 


.  in  ^  14,16,20,42-49 

effects  and  the  laser  pulse  can  be  optically  guided. 

index  of  refraction  associated  with  the  laser  in  the  plasma  is 

\  =  1  -  (l/2)(Wp/«)^  =  1  -  (l/2)(X/Xp)^Y"\/ngo,  (16) 

where  n^  =  n^^  +  Sn^,  Sn^  is  the  perturbed  electron  density,  and  y  is  the 
relativistic  factor  of  the  plasma  electrons.  For  a  long  laser  pulse  (long 
rise  time),  cT^  »  Xp,  it  may  be  shown  that  y  n^/n^^  =  (1  +  a^/2) 

A  necessary  requirement  for  optical  guiding  is  that  have  a  maximum  on 
axis,  an^/ar  <  O,  which  is  the  case  for  a  laser  intensity  profile  peaked  on 
axis.  Analysis  indicates  that  the  main  body  of  a  long  laser  pulse  will  be 
optically  guided,  provided  the  laser  power  P  exceeds  a  critical  value,  P  > 

P  ,  where^^'^^  P^[GW]  =  17(Xp/X)^.  As  an  example,  for  n^^  =10  cm  ,  Xp  = 
11  pm,  and  X  =  1  pm,  the  critical  laser  power  is  P^  =  2.1  TW.  For 
sufficiently  short  pulses,  <  l/«p’  plasma  has  insufficient  time  to 

respond  to  the  laser  pulse  and  relativistic  optical  guiding  does  not 
16,20 

occur. 

To  optically  guide  short  intense  laser  pulses  in  plasmas  a  preformed 
plasma  density  channel  can  be  used.l®'“  The  Index  of  refraction,  in  the 
absence  of  relativistic  effects,  is  given  by  Eq.  (16)  »ith  r  -  1  end  n^  . 
n  (r)  is  the  electron  density  radial  profile.  Optical  guiding  can  occur 
when  the  plasma  density  is  minimum  on-axis  and  analysis  of  the  wave 
equation  for  a  fixed  parabolic  plasma  density  channel  indicates  that 
optical  guiding  of  a  Gaussian  laser  pulse  occurs  when  the  channel  density 
depth  satisfies^®’^^  An^  >  An^f  where  An^  =  ^  critical 

channel  depth,  dn^  -  n^(r  -  r„)  -  n^(r  .  0)  and  r^  -  q'/mc^.  For  example, 
a  spot  size  of  r^  -  30  pm  implies  dn^  =  10*^  cm-^.  Simulations,  vhich 
include  the  self-consistent  evolution  of  the  density  channel,  indicate  that 
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density  channels  can  guide  short  pulse  lengths  and  high 
intensities. Recently,  the  guiding  of  modest  intensity  (10 
V/cmh  laser  pulses  in  a  plasma  density  channel  formed  by  an  axicon  focus 
has  been  experimentally  demons trated.^^  Capillary  discharges  can  also  be 
used  to  create  plasma  density  channels. Optical  beams  propagating  in 
plasmas  can  undergo  a  wide  range  of  instabilities  which  can  limit  the 
interaction  distance.  These  instabilities  include  backward  and  forward 
Raman  scattering, self-modulation, and  laser  hosing. 

C.  Self-Modulated  Laser  Wakefield  Acceleration 

In  the  recently  proposed  self-modulated  enhanced  electron 

acceleration  is  achieved  via  resonant  self-modulation  of  the  laser 

pulse. This  occurs  when  the  laser  pulse  length  is  several  plasma 

wavelengths,  and  the  peak  laser  power  satisfies  P  >  P^,  where 

P  [GW]  s  17(X  /X)^  is  the  critical  power.  At  fixed  laser  parameters,  both 

conditions  can  be  met  by  choosing  a  sufficiently  high  plasma  density,  since 

P  -  1/n  and  X  -  1/n^^^.  Wakefields  generated  in  the  self-modulated 
‘c  eo  p 

regime  are  more  than  an  order  of  magnitude  greater  than  those  generated  by 
a  laser  pulse  with  cx^  =  Xp/2,  assuming  fixed  laser  parameters. 

Acceleration  is  enhanced  for  four  reasons:  (i)  a  higher  density  produces  a 
larger  wakefield,  ~  n^^^,  (ii)  the  resonant  mechanism  excites  larger 
Wakefields  than  in  the  standard  LWFA,  (iii)  since  P  >  P^.,  relativistic 
focusing  further  enhances  the  laser  intensity,  and  (iv)  a  portion  of  the 
pulse  will  remain  guided  over  several  diffraction  lengths,  extending  the 

acceleration  distance. 

The  self-modulation  mechanism  can  be  understood  by  considering  a  long 
laser  pulse  cr^^  »  X^,  with  power  P  >  P^,  such  that  the  body  of  the  pulse 
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is  relativistically  guided.  The  leading  edge  of  the  pulse  will  create  a 

low-amplitude  wakefield  within  the  remainder  of  the  laser  pulse.  In  the 

Wakefield,  each  region  of  decreased  density  acts  as  a  local  plasma  channel 

to  enhance  the  relativistic  focusing  effect,  while  each  region  of  increased 

density  causes  defocusing.  This  results  in  a  low-amplitude  modulation  of 

the  laser  pulse  at  X„.  The  modulated  laser  pulse  resonantly  excites  the 

P  55 

wakefield  and  the  process  continues  in  an  unstable  manner.  This 

instability  develops  on  a  time-scale  associated  with  laser  diffraction.  In 

the  one-dimensional  limit,  laser  pulse  modulation  can  also  occur  via  the 

52  54 

usual  forward  Raman  scattering  (FRS)  instability. 

The  parameters  used  in  the  following  fluid  simulations  are  similar  to 
the  design  parameters  for  the  NRL  LWFA  experiment.  A  Gaussian  laser  pulse 
was  assumed  with  X  =  1  ym,  a^  =  0.96,  r^  =  10  ym  and  ct^^  =  120  ym  (400  fs), 
such  that  =  0.031  cm  and  P  =  2  TW.  The  density  is  n^^  =  1.1  x  10  cm 
(X  =  10  ym),  such  that  =  1.67  TW  (P  =  1.2P^).  At  t  =  0  the  laser  pulse 
is^outside  the  plasma  and  the  plasma  density  is  "ramped  up"  to  reach  full 
density  at  cT  =  0.1  cm  (3.2Zj^).  The  laser  pulse  is  initially  converging 
such  that  in  vacuum  it  would  focus  to  a  minimum  spot  size  of  r^  =  10  ym  at 

the  point  where  the  plasma  reaches  full  density  (0.1  cm).  The  simulation 

continues  until  ct  =  0.2  cm  (10  Zj^). 

Figure  5  shows  the  normalized  laser  intensity  \a^\  at  ct  =  0.2  cm  (at 
focus  in  vacuum,  =  aj).  Several  peaks  in  the  laser  intensity  are 

observable  with  period  ~  Xp.  The  corresponding  plot  of  the  axial  electric 

field  is  shown  in  Fig.  6.  As  the  laser  becomes  fully  modulated,  the 
resonant  excitation  of  the  field  becomes  stronger,  with  the  peak 
accelerating  field  increasing  Iron  |Ej  -  270  GV/m  at  cr  -  0.15  c.  to  IeJ 
=  350  GV/ra  at  ct  =  0.2  cm. 
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To  study  the  trapping  and  acceleration  of  beam  electrons  by  the 
Wakefield,  a  test  particle  code  is  used  to  accelerate  a  distribution  of 
non-interacting  particles  in  the  time-resolved  electric  and  magnetic 
Wakefields  of  the  fluid  simulation.  A  continuous  electron  beam  is  assumed 
with  initial  energy  of  1.0  MeV  and  normalized  emittance  =  7  mm-mrad. 

The  beam  is  initially  converging  such  that  in  vacuum  it  would  focus  to  a 
minimum  RMS  radius  rj,  =  350  ym  at  ct  =  0.1  cm.  After  ex  =  0.2  cm,  a  small 
fraction  of  the  original  particle  distribution  has  been  trapped  and 
accelerated.  Figure  7,  which  shows  the  peak  particle  energy  versus  ct, 
gives  a  peak  energy  of  113  MeV  at  ct  =  0.2  cm. 

V.  LASER-PLASMA  ACCELERATION  EXPERIMENTS 

Laser-plasma  acceleration  experiments  have  been  performed  in  a  number 

of  configurations  and/or  regimes.  For  example,  the  PBWA  uses  two  long 

pulse  (cTj^  »  X  )  laser  beams  of  frequencies  (o^  and  such  that  -  (<>2  = 

„  ,  to  resonantly  drive  a  plasma  .ave.‘2.38-«  The  LWFA,  on  the  other 
P  12-23 

hand,  uses  a  single,  intense  laser  pulse  to  drive  a  plasma  wave.  The 

LWFA  can  operate  in  the  short-pulse  (standard)  regime  in  which  CT^^  = 

X  ,  or  in  the  long-pulse  (self-modulated)  regime^^  in  which  cTj^  >  Xp. 

In  the  self-modulated  LWFA  regime,  laser  pulse  modulation  can  occur  by  the 

55  57 

FRS  instability^^ and/or  the  envelope  self-modulation  instability. 

Laser  acceleration  of  electrons  in  plasmas  has  been  observed  in 

59 

several  experiments  world-wide.  In  an  early  experiment  by  Joshi  et  al., 
a  single,  long  pulse  CO2  laser  beam  interacting  with  a  thin  Carbon  foil  was 
observed  to  produce  1.4  MeV  electrons  via  the  FRS  instability.  The 
pioneering  PBWA  experiments  at  UCLA^°  observed  a  28  MeV  energy  gain  (2.8 
GeV/m)  of  an  injected  electron  beam  (2  MeV)  using  two  CO2  laser  lines  in  a 
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plaSM  o£  density  9  x  These  experiments  are  partienlariy  eell- 

diagnosed  and  various  laser-plasaa  phenomena  have  been  observed.  PB»A 
experiments  at  Ecole  Polytechnique*^  observed  a  1.4  MeV  energy  gam  (0.6 
GeV/m)  o£  an  injected  electron  beam  (3.4  MeV)  using  t.o  Nd  laser  lines  in  a 
plasma  of  density  2  x  lo"  c.-’.  Experiments  at  Osaka=®  observed  electron 
acceleration  in  both  the  standard  and  self-modulated  LVFA  experiments  using 
a  8  TV,  1  ps  laser  pulse  and  a  source  of  -  1  MeV  electron  from  a  laser- 

solid  interaction.  In  the  standard  LVFA,  an  energy  gain  of  -  1  MeV  vas 

,  .  T  Joncitv  X  10^^  cm~^;  and  in  the  self-modulated 

observed  in  a  plasma  of  density  J  x  lu  cm  , 

LWFA,  an  energy  gain  of  17  MeV  was  observed  in  a  plasma  of  density  10 
cm-^  At  LLNL,^^  background  plasma  electrons  were  trapped  and  accelerated 
to  2  MeV  (~  2  GeV/m)  using  a  5  TW,  600  fs  laser  pulse  in  a  plasma  of 
density  2  x  10^^  cm'^  in  a  self-modulated  LWFA  experiment.  At 
Rutherford,^®  self-modulated  LWFA  experiments  generated  45  MeV  (~  90  GV/m) 
self-trapped  plasma  electrons  using  a  20  TW.  1  ps  laser  pulse  in  a  10^ 
cm-3  plasma.  Preliminary  results  from  Michigan^^  report  the  acceleration 
of  self-trapped  plasma  electrons  in  a  self-modulated  LWFA  experiment. 


VI.  DISCUSSION 

Several  important  issuea  aaaociated  .ith  laser  driver  acceleration  in 
vacuum,  neutral  gases  and  plasmas  have  been  discussed.  Laser  vacuum 
acceleration  is  limited  by  electron  slippage  and  diffraction  effects.  The 
slippage  distance  is  given  by  Eg.  (2).  A  general  class  of  laser  vacuum 
acceleration  configurations  requires  reflecting  surfaces  with  apertures  for 
the  electron  beam.  T»o  additional  limitations  apply  to  this  class  of 
configurations.  Material  damage  considerations  place  limits  on  the  laser 
intensity  electron  energy,  current  and  aperture  radius.  As  the  electron 
beam  passes  througb  the  aperture  a  portion  of  the  self-fields  is  reflected. 


20 


This  reflected  self-field  energy  represents  a  loss  in  electron  kinetic 
energy  and  places  a  limit  on  the  electron  energy  and  current.  Another 
general  class  of  laser  vacuum  accelerators,  e.g.,  the  vacuum  beat  wave 
accelerator,  utilizes  the  nonlinear  (ponderomotive)  laser  forces. 

The  electron  slippage  limitation  can  be  removed  by  using  a  neutral  gas 
for  the  accelerating  medium.  However,  diffractive  spreading  remains  an 
important  restriction.  To  overcome  both  diffraction  and  slippage,  a  self- 
guided  inverse  Cherenkov  acceleration  has  been  proposed  and  analyzed. 
Self-guiding  is  the  result  of  a  balance  between  diffraction,  nonlinear 
self-focusing  and  ionization  defocusing.  The  envelope  of  the  self-guided 
optical  beam  is  inherently  unstable  due  to  an  ionization-modulation 
instability.  Accelerating  gradients  of  -  0.5  GV/m  can  be  achieved  in  the 
self-guided  ICA. 

some  important  features  of  laser  acceleration  in  plasmas  have  been 
reviewed,  including  the  LWFA,  optical  guiding  in  plasmas,  and  the  self- 
modulated  LWFA.  A  simulation  of  the  self-modulated  LWFA,  based  on  the 
design  parameters  of  the  NRL  experiment,  indicates  that  energy  gains  in 
excess  of  100  MeV  and  gradients  in  excess  of  300  GV/m  are  possible. 

Several  laser-plasma  experiments  have  demonstrated  electron  acceleration, 
with  electron  energies  as  high  as  45  MeV  and  gradients  as  high  as  90  GV/m. 
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Fig.  1:  Schematic  of  electron  acceleration  in  vacuum  by  reflecting 
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Fig.  1:  Schematic  of  electron  acceleration  in  vacuum  by  reflecting 
higher  order  Gaussian  beam  from  a  mirror  placed  near  focus. 
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Fig.  2;  Acceleration  axial  field  -E^  plotted  versus  position  along  the  z- 
axis  for  the  higher  order  Gaussian  mode.  The  solid  curve  is  the 
high  energy  limit  y  Y  »  and  the  dashed  curve  is  the  low  energy 
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Laser  intensity  |  on-axis»  versus  C  at  ct  =  0.2  cm  in  the 
self-modulated  LWFA.  The  NRL  LWFA  design  parameters  were  used  in 
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Fig.  6:  Axial  electric  field  E  versus  C  at 
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CT  ^cm; 

Peak  electron  energy  versus  ex  for  a  beam  with  in 


